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We explore the scattering properties of ultracold ground-state Cs atoms at magnetic fields between 
450 G (45 mT) and 1000 G. We identify 17 new Feshbach resonances, including two very broad ones 
near 549 G and 787 G. We measure the binding energies of several different dimer states by magnetic 
field modulation spectroscopy. We use least-squares fitting to these experimental results, together 
with previous measurements at lower field, to determine a new 6-parameter model of the long-range 
interaction potential, designated M2012. Coupled-channels calculations using M2012 provide an 
accurate mapping between the s-wave scattering length and the magnetic field over the entire range 
of fields considered. This mapping is crucial for experiments that rely on precise tuning of the 
scattering length, such as those on Efimov physics. 



I. INTRODUCTION 

Cold cesium atoms have provided the foundation for 
many important experiments in basic science and also 
find application in precise atomic clocks. A thorough 
understanding of the collisions and interactions of two 
Cs atoms is crucial to interpret such experiments and 
optimize the applications. In particular, Cs has a com- 
plicated spectrum of magnetically tunable Feshbach res- 
onances [1], which allow precise control of the two-body 
interactions. These resonances are due to near-threshold 
bound states of the diatomic molecule Cs2 that can be 
tuned to match the near-zero energy of the colliding 
atoms. The resonances also allow an atomic sample to 
be converted with high efficiency into diatomic molecules 
by tuning an applied magnetic field across a resonance 
[2]. We have previously shown that the weakly bound 
molecules formed in this way can be converted into deeply 
bound [3] or even ground-state [4] molecules by stimu- 
lated Raman adiabatic passage. 

Early work with Cs atoms and its interactions was at 
low magnetic field, B < 250 G [5]. Results on low-field 
Feshbach resonances [6-9] made it possible to construct 
theoretical models of the near-threshold bound and scat- 
tering states of two cold Cs atoms [9, 10]. These models 
used the full Hamiltonian of Cs2 , including the potential 
energy curves of the singlet and '^S^ triplet states, 
the molecular spin-spin interaction, and the atomic hy- 
perfine interactions. Fitting the data allowed four key 
parameters of the model to be adjusted so that the res- 



onance structure could be reproduced accurately in the 
low-field region. These early models yielded an under- 
standing of the large clock shifts in Cs atomic fountain 
clocks [11], the anomalously large loss rates for colli- 
sions of doubly spin-polarized Cs atoms [12], and the 
magnetic field regions of moderate positive scattering 
length where Bose-Einstein condensation (BEC) was pos- 
sible [13]. Subsequent measurements of the binding en- 
ergies of weakly bound dimer states at fields up to 60 G 
[14, 15] were mostly in good agreement with calculations 
based on the model of Ref. [9] , which we designate M2004. 

Ultracold Cs is particularly important for the study 
of Efimov states [16], which are high-lying bound states 
of triatomic molecules that appear when the two-body 
interaction has a bound state very close to threshold. 
Efimov states cause additional loss features close to 2- 
body Feshbach resonances. The first observation of an 
Efimov resonance in Cs [17] was at a field near 8 G. We 
have recently extended this to observe Efimov features 
associated with three additional 2-body Feshbach reso- 
nances, at fields up to 900 G, and reached the remarkable 
conclusion that the Efimov features all occur at almost 
the same 2-body scattering length a, and thus all have 
almost the same binding energy at unitarity (a = oo) 
[18]. This universality of Efimov states was quite unex- 
pected, and contradicted previous theoretical indications 
[19], although subsequent theoretical work is now start- 
ing to explain its origins [20-23]. Efimov states in other 
ultracold systems are now also being found to show the 
same universal behavior [24-28]. 
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The new Efimov resonances used to demonstrate uni- 
versality in Ref. [18] are in the vicinity of two open- 
channel-dominated s-wave Feshbach resonances that 
were theoretically predicted with pole positions around 
550 G and 800 G [1]. Model M2004 is quite accurate 
at fields below 100 G, but at higher fields its predictions 
are in error by up to sa 1 per cent, or 8 G for B on the 
order of 800 G. This is not accurate enough to describe 
the scattering physics to sufficient precision to interpret 
the Efimov resonances, so an improved theoretical model 
is essential. 

The aim of the present paper is to describe the pre- 
viously unexplored properties of ^^'^Cs in its absolute 
atomic ground state in the magnetic high-field region be- 
tween 450 G and 1000 G. We present new measurements 
of both resonance positions and binding energies at fields 
up to 1000 G. We then use these results to develop a new 
optimized theoretical model, which we call M2012, that 
is accurate at both high and low magnetic fields. This 
model provides an accurate predictive tool to map the 
scattering length a{B) as a function of magnetic field B, 
which it is difficult to measure directly. This mapping 
was key to interpreting the three new Efimov features 
reported in Ref. [18]. 

This paper is organized as follows. Section II describes 
the essential molecular physics of threshold Cs states and 
describes the notation we use. Section III describes our 
experimental methods and results in the high-field region. 
Section IV describes our theoretical model. Section V de- 
scribes our least-squares fits to the experimental results, 
and compares experiment with theory in the regions of 
high, middle, and low fields. Section VI summarizes our 
conclusions. 



II. NEAR-THRESHOLD STATES OF CESIUM 
DIMER 

Two alkali-metal atoms in ^S states interact at short 
range to form singlet {X and triplet (a'^S^) states, 
with potential curves as shown in Figure 1. Levels that 
lie more than about 100 GHz below dissociation have 
fairly well-defined singlet or triplet character, so lie prin- 
cipally on one or the other of these curves. However, the 
levels of primary interest in the present work are very 
close to dissociation, and are bound by less than 1 GHz 
(and sometimes as little as 10 kHz). In this region the 
singlet and triplet states are strongly mixed by hyperfine 
interactions and it is more appropriate to describe the 
levels in terms of atomic quantum numbers. 

The zero-field levels of the Cs atom are characterized 
by the nuclear spin i = 7/2, the electron spin s = 1/2 
and their resultant / = 3 or 4, with the / = 4 level 
9.19 GHz above f — 3. In a magnetic field B, each level 
splits into 2/-I-1 sublevels labeled by m/, with the ground 
state {f,mf) = (3, -1-3). In the present work we label the 
atomic states with letters a,b,c, . . . in increasing order of 
energy. 
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FIG. 1. (Color online) Molecular potential energy curves 
V()(-R) and Vi(-R) for the singlet and triplet states of Cs2. 
The inset shows an expanded view of the long-range poten- 
tials separating to the two different / = 3 and 4 hyperfine 
states of the atoms at magnetic field _B = 0. 



For Cs2 there are 3 field-free atomic thresholds, la- 
beled in increasing order of energy by (/i,/2) — (3,3), 
(3,4) and (4,4), as shown in the inset of Fig. 1. In a mag- 
netic field, each threshold splits into sublevels labeled by 
(/ii JTii) -I- (/2, 7712). The near-threshold molecular states 
are to a good approximation described by quantum num- 
bers {fi,f2,F,Mp), where F is the resultant of /i and 
/2 and Mp = mi + 1712 (though mi and m2 are not 
individually conserved). Mp is a nearly good quantum 
number except near avoided crossings. For a homonu- 
clear molecule such as Cs2, F is also nearly conserved in 
the region where the atomic Zeeman effect is near-linear. 
Additional quantum numbers are needed for the molecu- 
lar vibration n and the partial-wave angular momentum 
L. For near-dissociation levels it is convenient to specify 
n with respect to dissociation, so that the topmost level 
is n = —1, the next is n = —2, etc. In the present work 
we describe near-threshold levels using a set of quantum 
numbers n{fif2)FL{Mp), with L = 0, 2, 4, etc. indicated 
by the usual labels s, d, g, etc. This is sometimes abbre- 
viated to FL{Mp) to avoid repetition. Following Ref. [1], 
we speak of bound levels with dominant s character in 
their wavefunction as s-wave levels; similarly for d- or g- 
wave levels with dominant L = 2 or 4 character in their 
wavefunctions. 

Each molecular level lies within a "bin" below its as- 
sociated threshold, with the boundaries of the bins de- 
termined by the long-range forces between the atoms. 
For Cs2, with V{r) = —Cer~^ at long range and Cg ~ 
6890 Ey^GQ for both the singlet and triplet states, the n — 
— 1 level lies between zero and —105 MHz, and the n = 
-2 level lies between -105 MHz and -725 MHz [1]. Sim- 
ilarly, bin boundaries can be worked out for more deeply 
bound levels. For Cs2 the background scattering length 
for each channel is large and positive, on the order of the 
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scattering length of the triplet potential. Under these cir- 
cumstances all the levels lie near the top of their respec- 
tive bins, and their energies E are approximately given 
by those of the triplet Born-Oppenheimer potential. Nu- 
merically, E/h is -0.0046 GHz, -0.11 GHz, -0.75 GHz, 
-2.4 GHz, -5.5 GHz, -10.6 GHz, and -18.1GHz for 
n -1 to -7 for the M2004 model [9]. 

Feshbach resonances occur where a weakly bound state 
exists at the same energy as the colliding atoms. Zero- 
energy Feshbach resonances thus occur at magnetic fields 
where a bound state crosses an atomic threshold. Each 
resonance is labeled by the quantum numbers of the 
bound state that causes it. We will work here with the 
scattering and bound states associated with the aa en- 
trance channel, with two atoms in state a with (/i, mi) — 
(3, -1-3). The energy zero at any magnetic field strength 
B is set to the energy of two a-state atoms. In s-wave 
scattering, the projection of the total angular momentum 
onto the field, Aftot, is thus always -1-6 and is a rigorously 
conserved quantity. The left-hand panels of Figure 2 
show the s-wave bound states with A/tot = Mp = -1-6 at 
magnetic fields up to 1000 G, together with the scattering 
length calculated using only s functions and thus includ- 
ing only resonances due to bound states with L — 0. 

The (/i,/2) = (3,3) levels with Mp = +6 in Fig. 2 
have the same magnetic moment as the separated aa 
atoms, and thus the energies of these bound states are 
parallel to the E — axis. However, Mp = +6 levels 
arising from other (/i,/2) combinations, and (3,3) levels 
with Mp 7^ -f6, have different magnetic moments and 
can cross the aa threshold as the magnetic field is varied. 
The three strong s-wave resonances at fields below 1000 
G are associated with ramping n{fif2)FL{Mp) states 
of-7(44)6s(6), -6(34)7s(6) and -6(34)6s(6) character. 
However, it should be noted that in each case the ramping 
s-wave state mixes strongly with the least-bound state 
— l(33)6s(6), which has a binding energy near 50 kHz, 
and this mixed state crosses threshold (and causes a pole 
in the s-wave scattering length) at a magnetic field below 
the field where the unperturbed ramping s-wave state 
would cross threshold. Such shifts in pole position are 
discussed in Refs. [1, 29, 30]. 

Figure 2 extends to negative magnetic field. This is 
to be interpreted as a reversal of axis, which is equiva- 
lent to changing the sign of all spin projection quantum 
numbers. Thus, the aa channel at —\B\ is equivalent to 
the gg channel at -I- 1 -B | , where the gg channel has two g- 
state atoms [(/i,toi) — (3,-3)]. The bound states and 
scattering length are continuous across -6 = 0, and in 
particular the low-field behavior of the scattering length 
is largely due to the ramping — 7(44)6s(6) state, which 
actually produces a resonance around B = —12 G (i.e., 
in the gg channel), as shown in detail in Fig. 3. 

Each s-wave bound state has a corresponding d-wave 
state, also with Mp = -1-6, that lies almost parallel to 
it but is shifted by the rotational energy of the vibra- 
tional state concerned; the rotational energy increases 
with binding energy and thus depends strongly on the vi- 



brational quantum number n. However, levels with L > 
and projection Ml such that Mtot = Mp + M^ = +Q will 
also cross the aa threshold and can contribute to s-wave 
threshold scattering [31]. The lower-right panel of Fig. 2 
shows the bound-state energies including the additional 
d-wave levels with Mp ^ 6. The upper-right panel shows 
the s-wave scattering length obtained with a basis set in- 
cluding both s and d functions (which we refer to as an 
sd basis set), showing the additional resonances that oc- 
cur. Figure 3 shows an expanded view of the scattering 
length and near-threshold bound states calculated with 
an sd basis set in the low-field region between —60 and 
-1-60 G. Figure 4 shows the near-threshold g-wave bound 
states, as studied in Refs. [9, 14, 31]. 

Previous work on ultracold Cs has focused on the 
low-field region. The first experimental studies of the 
coUisional properties of Cs were performed at Paris 
[12, 32, 33], Oxford [34, 35] and Stanford [6-8]. Chin 
et al. [7] and Leo et at. [10] studied over 30 resonances 
in several spin channels at fields below 130 G. Chin et 
al. [9] observed over 60 resonances in 8 different spin 
channels (aa, gg, hh, an, ao, ap, hn and gf) for fields 
up to 250 G. In recent years, we have explored the en- 
ergy spectrum of weakly bound Cs2 Feshbach molecules 
by magnetic moment, microwave [14] and magnetic field 
modulation spectroscopy [36]. The zero crossing of the 
s-wave scattering length was also precisely determined, 
using an approach based on measuring the interaction- 
induced dephasing of Bloch oscillations [37, 38]. 



III. FESHBACH SPECTROSCOPY AT HIGH 
MAGNETIC FIELD 

In the present work, we carry out a variety of differ- 
ent experiments on ultracold ^■^'^Cs in its lowest inter- 
nal state, (/, m/)=(3,3), at magnetic fields in the range 
between 450 and 1000 G. We first discuss the main ex- 
perimental procedures and conditions (Sec. HI A). Next 
we report on trap-loss spectroscopy, which allows mea- 
surement of the positions of narrow Feshbach resonances 
(Sec. HI B). Finally we present magnetic field modula- 
tion spectroscopy for the precise determination of molec- 
ular binding energies near broad Feshbach resonances 
(Sec. HI C). 



A. Sample preparation 

To access the high-field region, we have implemented a 
new magnetic-field system in the experimental setup [39]. 
This system is able to reach maximum magnetic field 
strengths up to 1400 G in a steady-state condition with 
10 mC long-term stability. The high magnetic bias fields 
are created by three separately controllable pairs of mag- 
netic field coils, made of 4 mm and 6 mm square-profile 
copper tubes insulated by glass-fiber braided sleevings. 
For each coil pair the electric currents, which are up to 
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FIG. 2. (Color online) Scattering length and energy levels versus magnetic field B for Cs collisions in the lowest energy aa spin 
channel. Panel (a): scattering length, calculated including only s basis functions. The vertical lines indicate the pole positions. 
Panel (b): Cs2 s-wave bound-state energies below the aa threshold with Mtot ~ +6, Mp — +6 in the range B = —200 to 
1000 G, calculated with s basis functions only. Panel (c): scattering length, calculated including s and d basis functions and 
including all matrix elements of the spin-dipolar coupling. Panel (d): Cs2 s-wave and d-wave bound state energies below the aa 
threshold with Mtot ~ +6 and all allowed values of Mp in the range B = —200 to 1000 G. The legend shows the L{Mf) labels. 
Levels with different L or Mp cross since small off-diagonal spin-dipolar matrix elements coupling them were not included in 
the calculation. Properties in the aa channel at negative values of B apply to the gg channel with reversed Mtot = —6 (see 
text). The black lines in the upper left corner of Panel (d) show the energies of the fg, ff, and eg atomic channels. 



400 A (4 mm tube) and 800 A (6mm tube), respectively, 
are supplied by two 6 kW power supplies, connected in 
parallel. The temperature of the coil system is kept be- 
low 50° C by internal water cooling of the copper tubes 
using a 10-bar pump system. Magnetic field stability is 
governed by controlling the current in the coils by an ac- 
tive feedback system, which operates at a precision level 
of 1Q~^. For this, the actual currents are measured by 
highly sensitive current transducers. We have checked 
that precise current control is sufficient for magnetic field 
control to the limit given above. Other influences, such as 
thermal expansion of the copper coils, play a minor role. 
A detailed description of the magnetic-field coil system 
can be found in Ref . [39] . 

The procedure used to prepare an ultracold cesium 



sample in the absolute atomic ground is based on well 
established cooling and trapping techniques, which are 
similar to the ones described in Ref. [40] down to the 
/iK-regime. After Zeeman slowing and cooling in a 
magneto-optical trap, the atoms are loaded into a three- 
dimensional optical lattice created by four laser beams, 
where Raman-sideband cooling [41-43] is performed for 
3.5 ms. During this stage, where a small magnetic field 
of several hundred mG is applied, the atoms are cooled 
and spin-polarized into the absolute ground state. After 
Raman-sideband cooling the ensemble size amounts to 
1.5x 10^ atoms at a temperature of about 1 /iK. Then, the 
atoms are transferred into a large- volume far-off-resonant 
dipole trap [44] generated by two crossed 100 W CO2 
laser beams, featuring a waist of about 600 /xm each. As 
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FIG. 3. (Color online) Panel (a): the real (red, solid) and 
imaginary (blue, dashed) parts of the complex scattering 
length (see Section IV), calculated with an sd basis set. For 
S < 0, two-body relaxation is possible to channels fg, ff, 
or eg. The inset shows an expanded view of the pole region 
of the resonance near — 12 G indicating a maximum variation 
in scattering length of « ±10^ ao due to resonance decay. 
Panel (b) shows the bound states below threshold, calculated 
as in Fig. 2 so that levels of different L or AIf cross instead 
of showing avoided crossings. The Mf labels in the negative 
B regions are shown with reversed sign. 



the optical trap is not strong enough to hold the atoms 
against gravity, an additional magnetic levitation field of 
31 G/cm is applied [13]. 

We use two different schemes, which we refer to as 
Scheme A and Scheme B. Scheme A is based on evapo- 
ration and detection at low magnetic fields, similar to our 
previous work [14] . This scheme could be implemented in 
a simple way, but ramping up to the high probe field and 
back down to the detection field involves crossing sev- 
eral Feshbach resonances, which causes additional losses 
and heating. In the course of performing the present ex- 
periments, we developed an improved approach (Scheme 
B) that allows imaging in the magnetic high-field region 
and optimization of evaporative cooling at higher fields. 
In the following, both schemes are described in detail. 
Figure 5 illustrates the generic timing sequence for both 
schemes. 
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FIG. 4. (Color online) Low-field near-threshold bound levels 
of s, d, and g symmetry. As for Figs. 2 and 3, the avoided 
crossings between states of different L or Mp are not cal- 
culated because spin-dipolar coupling is omitted in this cal- 
culation. The low-field g-wsve levels are of —2(33) charac- 
ter, whereas the level marked x2g{2) is of mixed n(/i/2) = 
-2(33), -6(34), and -7(44) character. 



In Scheme A, evaporation is performed for 2 s at con- 
stant depth of the CO2 laser trap in the magnetic low- 
field region. This stage of plain evaporation results in 
sa 5 X 10^ atoms at a temperature slightly below 1 /zK. 
Then, the CO2 laser trap is spatially overlapped with a 
crossed dipole trap created by a 1064 nm fiber laser, with 
waists of 40 /im and 250 /im. To continue evaporation, 
the tightly focussed 40 /im beam is ramped down from 
60 mW to 3.5 mW within 6.5 s, while the intensity of the 
250 ixm beam is fixed at 400 mW. During this procedure, 
both CO2 laser beams are switched off, finishing the sam- 
ple transfer. The s-wave scattering length is large and 
positive during the evaporation sequence and is adjusted 
for the final evaporation step to a ~ 200 aq, correspond- 



ing to a magnetic field strength of B^ 



21 G. Efficient 



evaporation conditions are encountered at this field be- 
cause of an Efimov-related three-body recombination loss 
minimum [17]. In this way, we end up with lO"" thermal 
atoms at a temperature of 70 nK in the magnetic low-field 
region. 

When ramping over the broad Feshbach resonances, 
sizeable effects of three-body recombination are unavoid- 
able, even when applying the fastest possible ramp 
speeds. This causes direct recombination losses and ad- 
ditional heating [45], which can cause subsequent evapo- 
ration losses in the measurement process. To avoid the 
latter effect, we recompress the trap by increasing the 
intensity of the 40 /im beam by about a factor 10 be- 
fore the magnetic field ramp is carried out. In the final 
step of the sample preparation, the levitation field is de- 
creased to 8 G/cm. The mean trap frequency is about 
a; = 27r X 46(5) Hz, and the final sample contains 10^ 
thermal atoms at a temperature of 120 nK in the mag- 
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netic low-field region. 

The measurements are performed by linearly ramping 
from Bevap within a ramp time tramp, i = 10 ms to the 
probe fields -Bprobc in the magnetic high-field region. As 
described above, crossing of the broad s-wave Feshbach 
resonances leads to considerable heating of the sample 
and additional particle loss. We estimate the tempera- 
ture at -Bprobc to be between 150 nK and 200 nK with 
Scheme A. 

To determine the particle number, we linearly de- 
crease the magnetic bias field to zero (i?imag — 0) within 
tramp, 2 = 10 ms and carry out resonant absorption imag- 
ing. The temperature is obtained in time-of-flight ex- 
pansion measurements after release from the trap. The 
magnetic field strength is determined from measurements 
of the (3, 3) — >■ (4, 4) microwave transition frequency by 
applying the Breit-Rabi formula [46] . 

In Scheme B, after loading the atoms from the Ra- 
man lattice to the levitated CO2 laser trap, the magnetic 
field is linearly ramped to the magnetic high-field region 
within 10 ms. For measurements performed below 800 G 
the ramp ends at 561 G (a = 1090 ao), whereas for mea- 
surements above 800 G it ends at 970 G (a = 1140 ao). 

At this stage, the three-body recombination losses that 
are encountered while crossing the broad s-wave Fesh- 
bach resonances are limited because of the low density 
and the relatively high temperature (T ~ 1 /iK) of the 
sample. The temperature dependence of losses follows 
from the unitarity limitation of three-body recombina- 
tion rates [47] . To compensate for the small change in the 
magnetic moment that is encountered during the ramp 
as a consequence of the quadratic contribution to the 
Zeeman effect, the magnetic levitation field is adjusted 
simultaneously. Then, 2 s of plain evaporation result in 
« 5 X 10^ atoms at a temperature of about 1 /iK (similar 
to Scheme A). Then, the CO2 laser trap is spatially over- 
lapped with the 1064 nm crossed dipole trap as described 
in Scheme A. 

In the crossed dipole trap, we perform forced evap- 
oration by decreasing the laser intensity of the AO 
beam from 60 mW to 3.5 mW within 15 s. During this 
step, the intensity of the 250 /im beam is changed only 
slightly, from 400 mW to 300 mW. Both CO2 lasers are 
switched off during the first 5 s of evaporation to achieve 
an efficient transfer to the 1064 nm trap. The magnetic 
bias field is adjusted during the evaporation sequence for 
optimized elastic scattering conditions. The last evap- 
oration step of this sequence ends at Bcvap = 558.7 G 
(a ~ 700 ao) or Bcvap = 894 G (a 300 ao), respectively. 
Note that close to 894 G an Efimov-related recombination 
minimum is present [18], which apparently facilitates ef- 
ficient evaporation. After recompression and reshaping, 
leading to a mean trap frequency of w = 27rx 26(3) Hz, we 
end up with a non-condensed sample of between 5 x 10** 
and 10^ atoms at a temperature of about 50 nK. 

For Scheme B, no broad s-wave Feshbach resonances 
are crossed in the final magnetic field ramps to reach 
Sprobc, and therefore no noteworthy heating effects and 




FIG. 5. Timing sequence for the magnetic field ramps for 
trap-loss spectroscopy and binding energy measurements. 
The magnetic field strength B is linearly ramped from the 
final evaporation field -Bcvap within the ramping time tramp, 1 
to the probe field -Bprobc- After an experimentally optimized 
hold time thoid ranging between 0.2 and 1 s, the magnetic 
field strength is (linearly) ramped to the imaging field -Bimag- 
For Scheme A, iJovap and -Bimag are in the magnetic low-field 
region- The magnified segment visualizes -Bprobo, which is 
constant in time for trap-loss spectroscopy (solid line) and 
sinusoidally modulated for the determination of binding en- 
ergies (dashed line)- 



particle losses are observed. The particle number is de- 
termined by high-field imaging in the vicinity of the 
zero crossing of the broad s-wave Feshbach resonances 
at -Bimag = 556.4 G and -Bjmag — 887.5 G. The magnetic 
field ramps to and from -Bprobc involve linear changes of 
the magnetic field with tramp. 1 = tramp, 2 = 10 ms, as for 
Scheme A. The procedures for magnetic field calibration 
and temperature determination are the same as the ones 
in Scheme A. 



B. Trap- loss spectroscopy 

Trap-loss spectroscopy is a well established method 
based on the enhancement of coUisional losses in the 
vicinity of a Feshbach resonance [1]. For atoms in the 
absolute atomic ground state, where inelastic binary col- 
lision are energetically forbidden, three-body recombi- 
nation is the lowest-order loss process. In this process 
three colliding atoms recombine to a molecule and a free 
atom. Typically, the kinetic energy released far exceeds 
the trapping potential, leading to loss of the three parti- 
cles involved. The general a^'-scaling of three-body re- 
combination rates [45, 48-50] leads to a maximum in 
losses at the magnetic field position -Bmaxj corresponding 
to the divergence of a at the Feshbach resonance pole, 
and a minimum at the position -Bmin, close to the zero 
crossing of the s-wave scattering length. This allows the 
observed losses to be directly related to the positions and 
widths of the Feshbach resonances. 

In this Section, we will first report on experiments 
characterizing the scattering properties in the vicinity 
of the broad s-wave Feshbach resonances, as shown in 
Fig. 2(a), by performing a broad magnetic field scan with 
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FIG. 6. Results of trap- loss spectroscopy in the magnetic field 
region between 450 G and 1000 G, performed with Scheme B. 
The enhancements in losses result from broad Feshbach res- 
onances centered at 494 G, 549 G and 787 G. The resonance 
positions, which are derived from our theoretical model (see 
Sec. IV), are indicated by dashed lines, labeled according to 
the quantum numbers FL{Mf) of the molecular states that 
cause the resonances. The measurements are performed with 
a hold time of 500 ms. Narrow Feshbach resonances are not 
visible in this scan because of the large step size of about 2 G. 
A remaining atom fraction of 1.0 corresponds to 8x 10* atoms. 



large step size. Then, we decrease the step size to per- 
form detailed scans to identify and characterize narrow 
Feshbach resonances, which originate from states with 
higher rotational angular momentum (L > 0). 

Trap-loss spectroscopy is performed by recording the 
remaining atom fraction after a hold time ihoid at the 
probe field Sprobo- In general, we cannot exclude ad- 
ditional losses encountered within iramp.i and iramp,2 



during the magnetic field ramps B, 



B 



cvap 



B 



probe 



and 



probe 



^imag- However, for the characterization of 



narrow Feshbach resonances only a small magnetic field 
region is investigated, where variations in the initial atom 
number are negligible. Furthermore, the measurements 
are performed with <hoid » iramp,i, iramp,2, strongly hm- 
iting the effect of finite ramp times. 

The broad scan of the magnetic high-field region cov- 
ers a range from 450 to 1000 G, as shown in Fig. 6. This 
scan clearly shows two broad loss features around 550 G 
and 800 G, which can be assigned to the two high-field 
s-wave Feshbach resonances, as discussed in Sec. II (see 
upper panel Fig. 2). These measurements demonstrate 
the large width of the Feshbach resonance near 800 G. 
Because of the unitarity limitation of three-body recom- 
bination losses [47], it is not possible to determine i?,nax 
accurately for the s-wave resonances by trap-loss spec- 
troscopy. 

In the region around 500 G, no s-wave Feshbach reso- 
nance is expected, but the theoretical model predicts a 
series of closely adjoining d-wave Feshbach resonances, 
as shown in Fig. 2. One of these has a width of about 
5G, producing the broad loss signal around 495 G seen 
in Figs. 6 and 7. 

We perform high-resolution scans by decreasing the 
step size of the magnetic field scans to a few mG. The 
results of these scans are displayed in Fig. 7. We observe 
15 narrow loss features, which can be assigned according 



to the theoretical model given in Sec. IV to Feshbach res- 
onances originating from the coupling of the free atoms 
to molecular states with rotational angular momentum 
i > 0. 

We observe three narrow resonances that cannot be 
attributed to s-, d- or g-wave molecular states in the 
present model. They are found at 461.62 G, 557.45 G 
and 562.17 G. Our model, however, predicts the existence 
of Feshbach resonances stemming from i-wave molecular 
states (L = 6) in the magnetic field regions where we 
observe these features. The calculations are not accurate 
enough to establish an unambiguous assignment, but the 
match between experiment and theory nevertheless pro- 
vides strong evidence that this is the first experimental 
observation of i-wave Feshbach resonances. These reso- 
nances are discussed further in Sec. IV. 

The positions of the poles of the d-, g- and i-wave 
Feshbach resonances obtained in these measurements are 
summarized in Table I. The peak positions are deter- 
mined by Gaussian fits to the loss peaks. For several of 
these resonances we also identify recombination loss min- 
ima, which also provide estimates for the corresponding 
resonance widths. 

Figure 8 shows expanded views of two regions in Fig. 7, 
where interesting cases of overlapping FR scenarios oc- 
cur. Figure 8(a) shows a (/-wave resonance centered at 
554.06(2) G sitting on the shoulder of the 549 G s-wave 
Feshbach resonance at a background scattering length of 
about — lOOOflQ. The zero crossing of a{B) leads to a 
loss minimum at 553.73(2) G. Even more intriguing is a 
d-wave Feshbach resonance situated at 820.4(2) G. There, 
the broad 787 G s-wave Feshbach resonance leads to an 
extremely large background scattering length of about 
— 4200ao. Experimentally, this large background masks 
the loss maximum but clearly reveals the zero crossing at 
819.41(2) G, as shown in Fig. 8(b). Both the g-wave and 
d-wave resonance have rather large widths, of 0.33(3) G 
and 0.96(22) G, respectively. Efimov-related three-body 
physics has been revealed in the vicinity of these reso- 
nances, as reported in Refs. [18, 51]. 



C. Binding energy measurements 

Binding energy measurements of weakly bound dimer 
states provide a powerful additional tool to extract in- 
formation on the cesium interaction potentials and scat- 
tering properties. In particular, for the s-wave Feshbach 
resonances the exact positions of the poles are obscured 
by strong loss across a broad magnetic field range, but 
can be extracted accurately from binding energies. In 
the present work, we measure the binding energies by 
magnetic field modulation spectroscopy, a method which 
was introduced in Ref. [52]. This method is based on 
a sinusoidal modulation of the magnetic bias field and 
allows the creation of dimers starting from an ultracold 
atom sample. This leads to an observable loss signal due 
to fast atom-dimer relaxation when the modulation fre- 
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FIG. 7. (Color online) Detailed results of trap-loss spectroscopy in the magnetic high-field region. We observe 15 Feshbach 
resonances, stemming from d-, g- and i-wave molecular states. Measurements indicated by (A)-symbols are obtained by 
Scheme A, whereas the (•)-symbols refer to data points acquired with Scheme B. The poles of the Feshbach resonances, 
omitting the s-wave resonances, are marked with an arrow. The loss features at 498.1 G, 553.3 G, 554.7 G, 818.9 G and 853.1 G, 
which are indicated by "E", are related to Efimov loss resonances as reported in Refs. [18, 51]. Note that the data in the 
intervals [728,729.5] G and [980,1000] G, which are measured for different thoid, are multiplied by scaling factors of 0.5 and 0.9, 
respectively, to reproduce the overall behavior shown in Fig. 6. A remaining atom fraction of 1.0 corresponds to 8 x 10* atoms. 
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FIG. 8. (Color online) Expanded view of Fig. 7, showing the 
Feshbach resonances arising from the g(3) state (a) and the 
— 6(34)6ci(6) state (b), which overlap with the broad s-wave 
resonances. The filled and empty arrows mark the poles and 
zero crossings of the scattering length, respectively. The re- 
maining atom fraction at the zero crossing in (b) is limited 
due to losses during the magnetic field ramps. 



quency matches the binding energy of the dimers plus 
the small relative kinetic energy of the colliding atoms. 
Since the modulation of the magnetic field is parallel to 
the magnetic bias field, only transitions between states 
with the same projection quantum number of the total 
angular momentum are observed. This procedure has 
been successfully applied in several experiments [36, 52- 
57] to determine atomic scattering properties. 

The atom samples for the binding energy measure- 
ments are prepared according to Scheme B, as described 
in Sec. Ill A. At iJprobo: a modulation signal is applied for 
a variable duration of thoid ~ 0.1 to Is, in a frequency 
range of 50 to 1600 kHz and an amplitude between 0.5 
and 3G (see Fig. 5). The amphtude and the duration 
of the pulse are experimentally adjusted for each bind- 
ing energy measurement to optimize the signal-to-noise 



ratio. The signal is generated by a programmable fre- 
quency generator and subsequently amplified by a com- 
mercial 25 W radio-frequency amplifier, which drives the 
current in a separate set of coils and thereby creates the 
modulation of the magnetic bias field. 

The measurements are usually performed by varying 
the modulation frequency at a fixed Bprobc- Another 
possibility, however, is to scan i?probo while the modu- 
lation frequency is kept constant. The advantage of the 
latter approach is that it is less sensitive to atom losses 
caused by technical imperfections, such as resonance phe- 
nomena in the electric circuit that drives the transitions. 
We checked that the two methods give consistent results 
in our measurements. Figure 9 shows sample loss signals 
derived in a frequency scan (a) and a magnetic field scan 
(b). 

We studied the binding energies E'b of the high-field 
s-wave states and of several d-, g- and i-wave states, as 
shown in Fig. 10. For the s-wave states with E\^/h < 
200 kHz, we observed asymmetric line shapes resulting 
from the finite temperature of the samples. We in- 
clude this effect in our fitting routine using the line- 
shape model of Ref. [58]. For s-wave states with Ei^/h > 
200 kHz and for dimer states with higher rotational an- 
gular momentum, the binding energy has a strong depen- 
dence on the magnetic field. In these cases, the magnetic 
field noise and the field gradient that is applied to levi- 
tate the atoms broaden and symmetrize the loss signals. 
For these symmetrized signals, the effect from the finite 
temperature plays a minor role, and we therefore obtain 
Eh by fitting a simple Gaussian distribution to the data. 

In the binding energy measurements, we observe sev- 
eral avoided crossings between molecular states. Around 
897 G and Ey,/h sa 500 kHz, the -6(34)6s(6) state crosses 
the — 6(34)5(i(5) state. These two states are clearly re- 
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molecular state 
n{hh)FL{Mp) 


J-Jmax \ ^ J 


J-Jmin l^*^ J 


-6(34)5d(5) 


460.86(5) 




i-wave'^' 


461.62(5) 




-6(34)7d(a;)2' 


492.45(3) 


492.63(3) 


-6(34)7d(a;)2' 


494.4(9) 


499.4(1) 


-6(34)7d(a;)2) 


501.24(3) 




-6(34)7d(a;)2' 


505.07(3) 






554.06(2) 


553.73(2) 


i-wave^' 


557.45(3) 




i-wave"^' 


562.17(3) 






565.48(3) 




-2(33)65(6)^) 


602.54(3) 






729.03(3) 




-6(34)6d(6) 


820.37(20) 


819.41(2) 


-6(34)5d(5) 


897.33(3) 




-6(34)6d(4) 


986.08(3) 





TABLE I. Results of trap-loss spectroscopy in the magnetic 
high-field region. The table shows the magnetic field values 
for loss maxima (Bmax), resulting from the poles of the Fesh- 
bach resonances, and minima (Bmin), which are related to the 
zero crossings of the s-wave scattering length. The numbers 
in brackets are the experimental uncertainties, including sta- 
tistical and systematic errors. The assignments of the Fesh- 
bach resonances identify the molecular states that cause the 
resonances. In the case of very narrow resonances, the zero 
crossings could not be determined experimentally. 
^'Por the i-wave resonances only the L-quantum number of 
the molecular state is known. 

^' These Feshbach resonances arise from d-wave molecular 
states with Mf = 4,5,6 and 7, which are strongly mixed 
at the atomic threshold. Therefore, we cannot give simple 
Mp quantum numbers and use "x" to indicate the strong 
coupling. 

^'For these states only the quantum numbers L and Mp are 
known. 

^^This molecular state is strongly mixed with the state 
-6(34)6(7(6). 



solved as separate loss features in each magnetic field 
scan performed at fixed frequency in the crossing region, 
as shown in the inset of Fig. 10(a). In addition, we ob- 
serve an avoided crossing at about 557 G and E\j/h « 
350kHz between the — 6(34)7s(6) state and an «-wave 
state, which clearly shows up in the binding energy mea- 
surements presented in Fig. 10(b). 

In view the large difference in the partial-wave angu- 
lar momentum AL = 6, the crossing around 557 G ap- 
pears to be surprisingly strongly avoided. To confirm 
this, we prepare molecular samples in the s-wave state 
by Feshbach ramps [1] and perform magnetic moment 
spectroscopy for magnetic field strengths ranging from 
556.5 G to 557.5 G using the Stern-Gerlach effect in the 
same way as described in Ref. [14]. To do this, we release 
the dimers from the trap while the magnetic field gradi- 
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FIG. 9. (Color online) Typical signals for magnetic field mod- 
ulation spectroscopy, (a) Frequency scan at a fixed magnetic 
field strength of 911.69 G. The asymmetric shape of the signal 
is fitted by a model that takes the temperature of the sam- 
ple into account [58]. (b) Magnetic field scan at a constant 
modulation frequency of 425 kHz, revealing a loss resonance 
at 557.02 G. This resonance stems from the s-wave molecular 
state. The line represents a Gaussian fit. A remaining atom 
fraction of 1.0 corresponds to (a) 4.5 x 10'' and (b) 2.5 x 10** 
atoms. 



ent is switched on. After a fixed time of flight, the dimers 
are dissociated by ramping back over either the 557.45 G 
or the 565.48 G Feshbach resonance. Subsequently, the 
atoms are imaged and the molecular magnetic moment 
is extracted from the vertical position of the atom cloud. 
We observe a smooth change of the magnetic moment 
around 557.15 G over a magnetic field range of about 
250 niG, indicating that the character of the molecular 
state also changes smoothly, from s-wave to i-wave char- 
acter, over the width of the crossing. Because of large 
calibration uncertainties in these measurements, we can- 
not provide absolute values for the molecular magnetic 
moments. In another experiment, we start with weakly 
bound s-wave Feshbach molecules at a magnetic field of 
560 G and attempt to jump the avoided crossing diabat- 
ically. As in the experiment previously outlined, we can 
simultaneously detect and distinguish s-wave and i-wave 
dimers by their magnetic moment. Applying a maximum 
ramp speed of 10 G/ms, the number of transferred dimers 
is below our detection limit. This sets a lower limit of 
30 kHz to the strength of the avoided crossing according 
to the Landau-Zener formula [59, 60]. 



IV. THEORETICAL MODEL 

The Hamiltonian for the interaction of two Cs atoms 
may be written 



2/z 



+ hi + h2 + V{R), (1) 



where /i is the reduced mass and is the operator for 
the end-over-end angular momentum of the two atoms 
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FIG. 10. (Color online) Results of experimental binding en- 
ergy measurements, (a) the — 6(34)6d(6) and — 6(34)5d(5) 
states (dashed lines) and the — 6(34)6s(6) state (solid line) 
between 820 G and 950 G. The lines are guides for the eye. 
The inset shows a binding energy measurement at the cross- 
ing of the -6(34)6s(6) and -6(34)5d(5) states at 509 kHz. At 
fields above 900 G the s-wave state takes on the character of 
the least-bound state —1(33) 6s (6), which has a binding energy 
around h x 50 kHz away from avoided crossings. The strong 
variation of the magnetic moment follows from the avoided 
crossing with this threshold s-wave state, (b) the binding en- 
ergies of the g{3) state (dashed line), an i-wave state (dot-dash 
line) and the — 6(34)7s(6) state (solid line). The positions of 
the resonance poles are marked by a (♦)-symbol, with values 
taken from Table I. 



about one another. The monomer Hamiltonians includ- 
ing Zeeman terms are 



(2) 



where ^ is the atomic hyperfine constant, si and §2 rep- 
resent the electron spins of the two atoms and ii and 12 
represent nuclear spins, and gn are the electron and 
nuclear g-factors, fiB is the Bohr magneton, and and 
iz represent the z-components of s and i along a space- 
fixed Z axis whose direction is defined by the external 
magnetic field B. The interaction between the two atoms 
V{R) is 



v{R) = y'=(i?) + y^(i?). 



(3) 



Here V'^iR) = Vo{R)V'-°'^ +Vi{R)V''^'^ is an isotropic po- 
tential operator that depends on the potential energy 
curves Vo{R) and Vi{R) for the respective X^E+ singlet 
and a'^E+ triplet states of Cs2, as shown in Figure 1. 
The singlet and triplet projectors V^°^ and V^^^ project 
onto subspaces with total electron spin quantum numbers 
and 1 respectively. The term V'^{R) represents small, 
anisotropic spin-dependent couplings that are responsible 
for the avoided crossings discussed in the present paper 
and are discussed further in Section IV C below. 



A. Computational methods for bound states and 
scattering 

The present work solves the Schrodinger equation for 
both scattering and Feshbach bound states of Cs2 by 
coupled-channel methods, using a basis set for the elec- 
tron and nuclear spins in a fully uncoupled representa- 
tion. 



\simsi)\iimii)\s2ms2)\i2mi2)\LML) 



(4) 



symmetrised to take account of exchange symmetry. The 
matrix elements of the different terms in the Hamiltonian 
in this basis set are given in the Appendix of Ref. [61]. 
The only rigorously conserved quantities are the parity, 
(—1)^, and the projection of the total angular momen- 
tum, Mtot = 'nT'si + "m-ii + + + Ml. The calcula- 
tions in this paper used basis sets with all possible values 
of TUs and for both atoms that are consistent with the 
required Mtot and parity, truncated at L,nax = 4 (an sdg 
basis set) unless otherwise indicated. All calculations in 
this paper are for s-wave incoming channels, so have even 
parity. 

Both scattering and bound-state calculations use prop- 
agation methods and do not rely on basis sets in the in- 
teratomic distance coordinate R. 

Scattering calculations are carried out using the 
MOLSCAT package [62] , as modified to handle collisions 
in magnetic fields [63]. At each magnetic field B, the 
wavefunction log-derivative matrix at collision energy E 
is propagated from i?niin = 6 ao to i?mid = 20 oq using 
the propagator of Manolopoulos [64] with a fixed step size 
of 0.002 oo, and from i?mid to -Rmax — 4,000 ao using 
the Airy propagator [65] with a variable step size con- 
trolled by the parameter TOLHI=10~^ [66]. Scattering 
boundary conditions [67] are applied at -Rmax to obtain 
the scattering S-matrix. The energy-dependent s-wave 
scattering length a{k) is then obtained from the diagonal 
S-matrix element in the incoming L = channel using 
the identity [68] 



a{k) 



1 

ik 



1-5, 



00 



1 + 5, 



00 



(5) 



where k^ — 2^E/h?. 

In the vicinity of a resonance at the lowest atomic 
threshold, the scattering length as a function of magnetic 
field (at small fixed k) follows the functional form 



a{B) - flbg [1 - A/(S - B„s)] 



(6) 



The resonance pole position B^^s may be associated with 
the 3-body loss maximum at a field -Bmax- For a nar- 
row resonance (where the background scattering length 
Cbg does not vary significantly across the resonance) , the 
width A is conveniently obtained from the difference be- 
tween the positions of the pole and zero in a{B). Ex- 
perimentally, this corresponds to the difference in field 
between the loss maximum at i^max and the loss mini- 
mum at Bmin- We have extended MOLSCAT to provide 
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an option to converge on poles and zeroes of a{B), instead 
of extracting them from a fit to a grid of points. 

Weakly bound levels for Feshbach molecules are ob- 
tained using a variant of the propagation method de- 
scribed in Ref. [61]. The log-derivative matrix is prop- 
agated outwards from i?inin to i?niid with a fixed step 
size of 0.002 uq and inwards from i?„iax to -Rmid with a 
variable step size, using the same propagators as for scat- 
tering calculations. i?mid = 25 qq and i?max = 4, 000 ao 
were used for most bound states, although i?mid = 35 ao 
and i?inax = 8, 000 Oq were needed for states within about 
50 kHz of dissociation. In Ref. [61], bound-state energies 
at a fixed value of the magnetic field B were located 
using the BOUND package [69], which converges on en- 
ergies where the smallest eigenvalue of the log-derivative 
matching determinant is zero [70] . However, for the pur- 
poses of the present work we used a new package, FIELD 
[71], which instead works at fixed binding energy and con- 
verges in a similar manner on the magnetic fields at which 
bound states exist. BOUND and FIELD both implement 
a node-count algorithm [70] which makes it straightfor- 
ward to ensure that all bound states that exist in a par- 
ticular range of energy or field are located. 

As described above, zero-energy Feshbach resonances 
can in principle be located as the fields at which the scat- 
tering length a{B) passes through a pole. However, with 
this method it is necessary first to search for poles, and it 
is quite easy to miss narrow resonances. However, since 
resonances occur at fields where there is a bound state at 
zero energy, the FIELD package provides a much cleaner 
approach: simply running FIELD at zero energy provides 
a complete list of all fields at which zero-energy Feshbach 
resonances exist. In the present work we located reso- 
nances using FIELD and then obtained their widths by 
converging on the nearby zero in a{B) using MOLSCAT. 



B. Representation of the potential curves 

At long range, the potentials are 

-C,oyw{R)/R'°±V,^{R), (7) 

where S = and 1 for singlet and triplet, respectively. 
The dispersion coefficients C„ are common to both po- 
tentials and the functions yn{R) account for retardation 
corrections [72]. The exchange contribution is [73] 



Fex(i?) = Aex(i?/ao)^exp(-^exi?/ao), 



(8) 



and makes an attractive contribution for the singlet and 
a repulsive contribution for the triplet. The value of /3cx 
is usually obtained from the ionization energies of the 
atoms [73], which for Cs gives /3cx — 1.069946, and 7 is 
related to /? by 7 = 7//3 — 1. In the present work we 
found that, to reproduce the experimental results, it was 
necessary to reduce /3ex slightly from its original value. 
We therefore introduce an additional factor p^x so that 
/3ox = 1.069946pcx, with 7 adjusted accordingly. 



The detailed shapes of the short-range singlet and 
triplet potentials are relatively unimportant for the ul- 
tracold scattering properties and near-threshold binding 
energies considered here, although it is crucial to be able 
to vary the volume of the potential wells to allow adjust- 
ment of the singlet and triplet scattering lengths. In the 
present work we retained the functional form used by Leo 
et al. [74] and Chin et al. [75]. Each short-range poten- 
tial is represented by a set of 14 ab initio points between 
R = 7 and 20 oq [76]. The two sets of potential points 
are first multiplied by R^ and the resulting (smoother) 
functions are interpolated using Akima splines [77] to ob- 
tain their values at i?LR — 17.6 uq. The value of A^x is 
chosen to match Vc^ to (Vi — Vb)/2 at Rlr, and both 
sets of points are shifted to match {Vi + Vo)/2 at i?LR,- 
Finally, the analytic Vg^{R) is used to generate new grid 
points between 17.6 ao and 20 oq and the resulting sets 
of points are reinterpolated as above between R = 7 and 
20 ao- The analytic long-range form (7) is used outside 
20 ao. 

The flexibility needed to adjust the singlet and 
triplet scattering lengths is provided by simply adding 
a quadratic shift to each of the singlet and triplet poten- 
tials inside its minimum. 



Vt'*\R)^SsiR 



Res)^ 



for 



R<ReS, (9) 



with i?eo = 8.75 ao and i?ei — 11.8 ao. 



C. Magnetic dipole interaction and second-order 
spin-orbit coupling 

At long range, the coupling V'^{R) of Eq. (3) has 
a simple magnetic dipole-dipole form that varies as 
1/R^ [78, 79]. However, for atoms as heavy as Cs, second- 
order spin-orbit coupling provides an additional contri- 
bution that has the same tensor form as the dipole-dipole 
term and dominates at short range [80]. In the present 
work, V'^{R) is represented as 

V^'iR) = A(i?) (Ji • S2 - 3(si • eR)i.S2 ■ cr)) , (10) 

where br is a unit vector along the internuclear axis and 
A is an i?-dependent coupling constant. The second-order 
term has been calculated by Kotochigova et al. [81] and 
fitted to a biexponential form, so that the overall form of 
A(i?) is 



X{R) = EW 



Afg-*exp (-2/32Soi?) 
-A'°^§exp{-(32soR)^ 



(i?/«o)= 



,(11) 



where a « 1/137 is the fine-structure constant and 
the parameters obtained from fitting to the electronic 
structure calculations [81] are A'^°Q^/hc — 34.4 cm~^, 
A^^gQ/hc = 0.25 cm~^ and /32so = —0.35 a^^. How- 
ever, in fitting to the experimental results, this coupling 
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function was found to be slightly too strong. We there- 
fore retained the functional form (11) but introduced an 
additional scaling factor S2S0 that multiplies both expo- 
nential terms and is allowed to vary in the least-squares 
fit to the experimental results. 



V. LEAST-SQUARES FITTING OF POTENTIAL 
PARAMETERS 

In the present work, our primary objective is to obtain 
potential parameters that give a reliable representation 
of a{B) in the regions where Efimov resonances occur, 
namely near 8 G, 554 G and 853 G. Earlier potentials 
[74, 75] focused on representing the positions of Feshbach 
resonances in the low-field region below about 60 G. 

A key advantage of the propagator approach to lo- 
cating bound states and resonances, implemented in the 
BOUND and FIELD programs, is that it is fast enough 
to be incorporated in a least-squares fitting program. We 
have therefore carried out direct least-squares refinement 
of the potential parameters. We experimented with fit- 
ting various combinations of parameters, and concluded 
that adequate flexibility is available in the 6-parameter 
space Sq, Si, Cg, Cs, £'230, Pex- 

The set of experimental results used for fltting is listed 
in Table II. It consists of all the observed resonance 
positions (loss maxima) in the regions relevant to Efi- 
mov physics, together representative bound-state posi- 
tions. It also includes zero crossings of the scattering 
lengths measured from loss minima, Bloch oscillations 
[37] and BEG collapse measurements [82]. All the bound- 
state positions are expressed as the fields at which bound 
states exist at specific binding energies, except for the 
— 3(33)6g(6) state, which is almost parallel to the lowest 
threshold and is included in the fit as a binding energy 
a.t B ~ 18.6 G. The strength of the avoided crossing be- 
tween the -2(33)4d(4) and -7(44)6s(6) states near 48 G 
is included as an explicit difference between fields where 
states exist at energies of 78 kHz and 24 kHz. The quan- 
tity optimized in the least-squares fits was the sum of 
squares of residuals ((obs— calc)/uncertainty), with the 
uncertainties listed in Table II [83]. 

Although the experimental data do allow all 6 of the 
potential parameters described above to be determined, 
the fit is very highly correlated. Under these circum- 
stances, a fully automated approach to fitting is unre- 
liable: individual least-squares steps often reach points 
in parameter space where the levels have moved too far 
to be identified reliably, particularly in the early stages 
of fitting. We therefore carried out the fitting using 
the I-NoLLS package [84] (Interactive Non-Linear Least- 
Squares), which gives the user interactive control over 
step lengths and assignments as the fit proceeds. This 
allowed us to converge on a minimum in the sum of 
weighted squares, with the parameters given in Table HI. 

The parameter uncertainties are given in Table HI as 
95% confidence limits [85]. However, it should be em- 



phasized that these are statistical uncertainties within 
the particular parameter set. They do not include any 
errors due to the choice of the potential functions. Such 
model errors are far harder to estimate, except by per- 
forming a large number of fits with different potential 
models, which is not possible in the present case. 

In a correlated fit, the statistical uncertainty in a fitted 
parameter depends on the degree of correlation. How- 
ever, to reproduce the results from a set of parame- 
ters, it is often necessary to specify many more digits 
than implied by the uncertainty. A guide to the number 
of digits required is given by the parameter sensitivity 
[85] , which essentially measures how fast the observables 
change when one parameter is varied with all others held 
fixed. This quantity is included in Table HI. 

The singlet and triplet scattering lengths and the pole 
positions of the s-wave resonances are not directly ob- 
served quantities. Nevertheless, their values may be ex- 
tracted from the final potential. In addition, the statisti- 
cal uncertainties in derived parameters such as these may 
be obtained as described in Ref. 85. The values and 95% 
uncertainties obtained in this way are give in Table IV. 

A. Region between 800 G and 920 G 

In the region between 800 and 920 G, the near- 
threshold molecular structure and the corresponding 
scattering properties are relatively straightforward. As 
seen in Fig. 2, the ramping — 6(34)6s(6) state has a strong 
avoided crossing with — 2(33)6s(6) near 850 G at a bind- 
ing energy of about 110 MHz. This crossing is still in- 
complete when the mixed state crosses the least-bound 
state — l(33)6s(6) around 100 kHz. It is the state result- 
ing from this second crossing that is observed at binding 
energies of 0.1 to 1.7 MHz between 890 and 920 G. These 
points are very well reproduced by the fit, as seen in Fig- 
ure 11. The — 3(34)6(i(6) state that crosses the axis at 
820.37 G is also well reproduced by the fit, as seen in 
Figure 12. 

The position of the s-wave resonance pole in this region 
has not been observed because of large 3-body losses, but 
it is predicted by the fitted potential at 786.8(6) G. It is 
interesting to note that, because of the shift produced 
by the crossing with the least-bound state, the s-wave 
resonance occurs at a lower field than the corresponding 
d-wave resonance in this case, even though (as always) 
the unperturbed ramping s-wave state is lower in energy. 

B. Region around 550 G 

The region around 550 G is considerably more com- 
plicated. In this case the crossing between the ramping 
-6(34)7s(6) state and -2(33)6s(6) near 510 G is much 
narrower, so is almost complete by the time — 6(34)7s(6) 
crosses the least-bound state —1(33) 6s (6). This produces 
a zero crossing near 556 G. However, as described in the 
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TABLE II. Quality of fit between calculations using the M2012 model and the experimental results used in the fit. 





-Bobs (G) Scale (G) -Bobs 


-^calc 


(G) Unc 


• (G) 


method 


reference 


-7(44)6s(6) at 7.8 MHz 


17 


53 


17 


51 


0.02 





02 


microwave spectr. 


[14] 


-7(44)6s(6) at 1.2 MHz 


21 


60 


21 


59 


0.01 





02 


microwave spectr. 


[14] 


-7(44)6s(6) at 104 kHz 


32 


05 


31 


70 


0.35 





03 


magnetic field mod. 


[36] 


Zero crossing near 17 G 


17 


12 


17 


14 


-0.02 





01 


Bloch oscillations 


[37] 


-2(33)4d(4) at 174 kHz 


48 


01 


48 


01 


0.00 





06 


magnetic field mod. 


[36] 


-2(33)4d(4) crossing strength 78-24 kHz 


1 


19 


1 


21 


-0.02 





02 


magnetic field mod. 


[31] 


Loss minimum (d) near 48 G 


47 


94 


47 


98 


-0.04 





04 


inferred from magn. field mod. 


[36] 


Loss maximum (d) near 48 G 


47 


78 


47 


79 


-0.01 





06 


inferred from magn. field mod. 


[36] 


A (d) near 48 G 





16 





18 


-0.02 





06 


inferred from magn. field mod. 


[36] 


2g{2) at 17 kHz 


53 


42 


53 


76 


-0.34 





08 


magnetic field mod. 


[36] 


-2(33)6p(6) at 18.6 G (MHz, not G) 


-5 


03 


-4 


99 


-0.04 





01 


microwave spectr. 


[14] 


Loss maximum — 6(34)7d(a;) 


492 


45 


492 


68 


-0.23 





06 


trap loss spectr. 


this work 


Loss maximum — 6(34)7d(a;) 


501 


24 


501 


44 


-0.20 





06 


trap loss spectr. 


this work 


Loss maximum — 6(34)7d(a;) 


505 


07 


505 


37 


-0.30 





06 


trap loss spectr. 


this work 


-6(34)7s(6) at 1.0 MHz 


556 


47 


556 


48 


-0.01 





02 


magnetic field mod. 


this work 


-6(34)7s(6) at 700 kHz 


556 


72 


556 


76 


-0.04 





02 


magnetic field mod. 


this work 


-6(34)7s(6) at 170 kHz 


557 


88 


557 


80 


0.08 





03 


magnetic field mod. 


this work 


-6(34)7s(6) at 100 kHz 


558 


44 


558 


36 


0.08 





03 


magnetic field mod. 


this work 


Zero crossing near 556 G 


556 


26 


556 


19 


0.07 





03 


collapse of BEG 


[82] 


-3(33)6p(3) at 368 kHz 


553 


44 


553 


44 


-0.00 





01 


magnetic field mod. 


this work 


Loss minimum (g) near 554 G 


553 


73 


553 


75 


-0.02 





01 


trap loss spectr. 


this work 


Loss maximum (g) near 554 G 


554 


06 


554 


07 


-0.01 





02 


trap loss spectr. 


this work 


A {g) near 554 G 





33 





32 


0.01 





01 


trap loss spectr. 


this work 


-6(34)6s(6) at 1.7 MHz 


890 


52 


890 


61 


-0.09 





02 


magnetic field mod. 


this work 


-6(34)6s(6) at 356 kHz 


899 


93 


900 


19 


-0.26 





03 


magnetic field mod. 


this work 


-6(34)6s(6) at 110 kHz 


915 


66 


915 


54 


0.12 





03 


magnetic field mod. 


this work 


Zero crossing near 881 G 


880 


90 


880 


66 


0.24 





03 


collapse of BEG 


[82] 


-6(34)6d(6) at 342 kHz 


819 


17 


819 


20 


-0.03 





03 


magnetic field mod. 


this work 


Loss minimum (d) near 820 G 


819 


41 


819 


37 


0.04 





03 


trap loss spectr. 


this work 


Loss maximum (d) near 820 G 


820 


37 


820 


33 


0.04 





02 


trap loss spectr. 


this work 


A (d) near 820 G 





96 





97 


-0.01 





05 


trap loss spectr. 


this work 



TABLE HI. Parameters of the fitted potential. 





fitted value 


confidence 


sensitivity 






limit (95%) 




Cq (-Ehao) 


6890.4768 


0.081 


0.0003 


Cg (-Ehflo) 


1009289.6 


2900 


0.1 


5*0 (-EhiQ ^) 


3.172749 X 10"* 


2.5 X lO"'' 


2.7 X 10"^" 




1.343217 X 10"'' 


2.5 X 10"" 


1.3 X 10"*" 


5*230 


1.36432 


0.022 


0.00065 


Pcx 


0.978845 


0.0026 


5.5 X 10"'' 



TABLE IV. Comparison between key quantities calculated 
from different potentials. 



derived 


M2012 


M2004 


parameters 






as (bohr) 


286.5(1) 


280.37(6) 


At (bohr) 


2858(19) 


2440(24) 


s pole near —10 (G) 


-12.38(8) 




s pole near 550 (G) 


548.78(9) 




s pole near 800 (G) 


786.8(6) 





experimental section, there is also a ramping —3(33)3(7(3) 
state that crosses threshold just below this, producing an 
additional pole and zero crossing near 554 G. Because of 
the large background scattering length arising from its 
proximity to the s-wave resonance, the (7-wave resonance 
is far wider than is usual. 



The bound-state measurements near 557 G are further 
complicated by an additional state that crosses and ap- 
pears to mix with the s-wave state in the top 1 MHz. 
This cannot be assigned as s-, d- or g-wave; it must be 
due to an i-wave state, and indeed the M2004 model 
predicts four i-wave levels to be in the range between 
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FIG. 11. (Color online) (a) The scattering length in the region 
between 875 G and 950 G. The solid lines show results calcu- 
lated for the M2012 potential with an sdg basis set. The blue 
diamond near 880 G indicates the position of the zero cross- 
ing, (b) The calculated bound-state energies with the same 
basis set. The blue diamonds show the measured energies in 
this region. 
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FIG. 13. (Color online) (a) The scattering length between 546 
G and 564 G. The solid lines show results calculated for the 
M2012 potential with an sdg basis set. The blue diamonds 
show the measured pole position and loss minimum near 554 
G and the zero crossing near 556 G. The dashed line shows the 
calculated pole position, (b) The calculated binding energies 
with the same basis set as for Panel (a). The blue diamonds 
show the measured binding energies. 
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FIG. 12. (Color online) The scattering length and bound- 
state energy in the pole region of the — 6(34)6ci(6) resonance. 
The solid lines show results from model M2012, calculated 
with an sdg basis set. The blue points in panel (a) show the 
measured pole position and loss minimum. The blue points 
in panel (b) show the measured binding energies. The blue 
points in panel (b) show the measured binding energies. The 
dashed line in (a) and the arrow in (b) shows the calculated 
pole position. 



around 530 G and 590 G, with two between 555 G 
and 565 G. However, these states are associated with 
hyperfine-excited thresholds, and their exact positions 
are very sensitive to details of the potential that do not 
significantly affect the other experimentally measured 
quantities considered here. We have been unable to de- 



cide unambiguously which one is responsible for the ob- 
served crossing. In addition, the j-wave states have very 
little influence on the s-wave scattering length that is 
the main object of interest in this region. We therefore 
decided to fit using a basis set with Lmax = 4, which 
excludes the z-wave states entirely, and also to exclude 
from the fit any binding energies affected by the crossing 
between the s- and i-wave states. 

The general fit to the ramping s-wave bound state is 
shown in Fig. 13 and an expanded view of the fir to the 
f;-wave state is shown in Fig. 14. It may be seen that the 
calculated s-wave state passes well through the experi- 
mental points either side of the crossing with the i-wave 
state, while the calculated g-wave state reproduces the 
bound-state energies as well as the zero crossing and pole 
in the scattering length. 

The crossing between the i-wave state and the 
— 6(34)7s(6) state near 557 G is surprisingly strongly 
avoided. As described above, we were unable to ramp 
the magnetic field fast enough to transfer a detectable 
number of molecules between the two states, which sets 
a lower limit of 30 kHz to the strength of the avoided 
crossing. In a zeroth-order picture, the two states are 
separated by AL — 6 and the only coupling off-diagonal 
in L is V*^ {R) , which can couple only AL = 2 and is quite 
small. However, quantitative modeling of this effect re- 
quires a theoretical model that places the i-wave state in 
the correct place. This is likely to require simultaneous 
modeling of the present results and the spectroscopy of 
more deeply bound levels, and this remains a subject for 
future work. 
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FIG. 14. (Color online) Expanded view of the scattering 
length in the pole region of the —3(33)6(7(3) resonance. The 
solid line shows results calculated for the M2012 potential 
with an sdg basis set. The points show the measured bound- 
state energies. The calculated (red) and measured (blue) pole 
positions are indicated by arrows. 



C. Low- field region 

The low-field region, below 60 G, is also quite com- 
plicated. The ramping s-wave state responsible for the 
Efimov resonances in this region is — 7(44)6s(6). How- 
ever, there are also families of ramping —2(33)6(7 states 
that cross threshold between 4 and 8 G and — 2(33)4(; 
states that cross between 11 and 21 G. In addition, there 
is a — 2(33)4(i(4) state that crosses near 48 G and a 2^(2) 
state that crosses near 54 G. The latter state does not 
carry a clear n(/i/2) signature, but it clearly has 2g{2) 
character. 

The binding energies of many of these states have been 
measured by magnetic moment spectroscopy at binding 
energies up to about 10 MHz [86], though these mea- 
surements have significant uncertainties associated with 
the integration over field. However, there are much more 
precise results for the — 7(44)6s(6) state as it crosses with 
the least-bound state, obtained from microwave spec- 
troscopy [14] and magnetic field modulation spectroscopy 
[36]. Particularly important are measurements of the 
crossing between this mixed state and the — 2(33)4(i(4) 
state near 48 G, since the strength of this crossing pro- 
vides the most direct experimental information available 
on the strength of the second-order spin-orbit coupling 
and thus on the potential parameter S'2SO- 

Figures 15 and 16 show the overall fit to the bound 
states below 60 G for both the M2004 and M2012 po- 
tentials. The ramping s-wave state between 17 G and 
60 G is of particular interest. This state switches from 
— 7(44)6s(6) to — l(33)6s(6) character as B increases. All 
the experimental energies are well reproduced, including 
those near the two avoided crossings with the 4(i(4) and 



25(2) states shown in Fig. 15. Even in this region, the 
M2012 model agrees with the experimental results signif- 
icantly better than the M2004 model. 
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FIG. 15. (Color online) The bound-state energies between 
35 G and 60 G. The solid lines show results calculated for 
the M2012 potential with an sdg basis set, showing one side 
of the avoided crossings of the — 7(44)6s(6) bound state with 
the — l(33)6s(6) least-bound state in the threshold entrance 
channel. The points show the previously measured bound- 
state energies [14, 36]. The agreement is much better with 
the new M2012 model than with the old M2004 model, shown 
as light dashed lines. 



D. Independent tests of the M2012 model 

The older potential models [74, 75] satisfactorily re- 
produced the bound states for s-, d- and g-wa.ve states 
at fields below 60 G. However, the measurements also 
revealed the existence of /-wave (L = 8) states [14, 15], 
which do not lead to observable Feshbach resonances be- 
cause of their very weak coupling to the s-wave thresh- 
old channel. The M2004 model failed for these /-wave 
states, with errors of up to 5 G in the positions at a 
given bound-state energy. The /-wave states were not in- 
cluded explicitly in our fits, but the comparison between 
the calculated levels and experiment is shown in Fig. 16 
for both the current model M2012 and the older M2004 
one. It is clear that M2012 gives a far more satisfactory 
reproduction of the experimental /-wave levels. 

There is a particularly interesting region near 500 G, 
where a group of four strongly coupled —6{3A)7d{Mp) 
levels cross threshold. The underlying levels have Mp 
values of 4, 5, 6, and 7, but are strongly mixed with 
one another so that Mp is not a good quantum number 
for the actual eigenstates. The bound states have not yet 
been measured in this region, but the calculated levels are 
shown in Fig. 17, together with the calculated scattering 
length and the positions of the measured loss maxima 
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FIG. 16. (Color online) Comparison between experimental 
results for s-wave levels (green diamonds [14, 36]) and g- 
and /-wave levels (open circles [14, 15]) with the results of 
the M2012 model (solid lines). Only levels for which experi- 
mental results exist are shown. The s-wave and g-wave levels 
are calculated with the full sdg basis set. The Z-wave levels 
are calculated with an I basis only. The light dashed lines 
show the Z-wave levels calculated with the M2004 model. The 
arrows show the separation between equivalent levels in the 
M2004 and M2012 models. The additional quantum numbers 
n(/i/2) are —2(33) for the two lig levels and —3(33) for the 
three 6Z levels. 



FIG. 17. (Color online) The scattering length (a) and bound- 
state energies (b) between 490 G and 508 G, calculated for the 
M2012 potential with an sdg basis set. The solid diamonds 
in (a) show the positions of experimental loss maxima and 
the open circles show the positions of loss minima. The four 
ramping — 6(34)7d(a;) levels, labeled collectively by x since 
their Mp components are mixed, undergo avoided crossings 
with the — 2(33)6(;(6) level that passes through this region 
with —E/h « 2 MHz. The avoided crossings with the least- 
bound level — l(33)6s(6) near —E/h « 0.01 MHz are too close 
to the E = line to be seen on the figure. 



and minima. The loss maxima are well reproduced. The 
two loss minima at 499.6 G and 502.15 G are also well 
reproduced by the model, but the loss minimum at 492.8 
G is not near a zero-crossing of the calculated scattering 
length. However, there are two strong overlapping and 
interfering resonances with poles in a{B) at 492.7 G and 
495.0 G, and it is not clear how to interpret the three- 
body loss in such a region. This complex region from 490 
G to 510 G needs further investigation, especially since 
it may display rich Efimov physics [51]. 



E. Mapping between scattering length and 
magnetic field 

An important goal of this paper is to develop a theo- 
retical model that is capable of giving an accurate map- 
ping between the scattering length a{B) and the mag- 
netic field with particular focus on collisions between 
two Gs atoms in the lowest Zeeman level of the ground- 
state manifold. Our new M2012 model is based on real- 
istic potentials and includes the full spin Haniiltonian for 
the Cs2 molecule at long range. It has been calibrated 
against experimental binding energies and a large num- 
ber of resonance positions at magnetic fields between 
and 1000 G. We therefore expect it to provide an accurate 
representation of a{B) across this entire range of fields. 
Figure 18 shows a{B) on a grid with 0.1 G spacing over 
the new experimental range between 460 G and 1000 G. 



At this resolution, some narrow resonances are not fully 
resolved, and the narrowest ones are not visible at all. 

The Supplemental Material provides a tabulation of 
a{B) against B over the full range of fields from to 
1200 G, for use in interpreting experiments. Additional 
tables lists all s-, d-, g-, i- and Z-wave resonances up to 
1000 G. 



VI. CONCLUSION 

We have explored the scattering properties of ground- 
state Cs atoms and the binding energies of weakly bound 
Cs2 molecules in the previously uncharted magnetic field 
range up to 1000 G, using a combination of experiment 
and theory. We have developed a new model of the in- 
teraction potential that reproduces the experiments ac- 
curately over the entire range of magnetic fields studied. 

Experimentally, we have investigated Feshbach reso- 
nances and dimer binding energies in the magnetic field 
range between 450 G and 1000 G, utilizing an ultracold 
Gs sample in an optical dipole trap. Around 550 G and 
800 G, we verified that the general scattering properties 
of atomic Gs are governed by two broad s-wave Feshbach 
resonances. Fifteen new Feshbach resonances stemming 
from molecular states with L > were pinpointed by 
trap loss spectroscopy. We found the first evidence for 
the existence of i-wave Feshbach resonances, resulting 
from the coupling of molecular states with L = 6 to the 
atomic threshold. By performing magnetic field modula- 
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FIG. 18. (Color online) The scattering length between 460 
G and 1000 G, calculated for the M2012 potential with an 
sdg basis set. The dashed lines indicate the pole positions of 
the broad s-wave resonances near 549 G and 787 G and the 
d-wave one near 820 G. Note that some narrow resonances 
appear as spikes because of the limited grid spacing of 0.1 G, 
and some of the narrowest are not visible at all. 

tion spectroscopy, we determined the binding energies of 
several dimer states, paying particular attention to the 
two s-wave states that are responsible for the general Cs 
scattering properties in the high-field region. 

To calculate the scattering properties and bound- 
state energies, we solve the Schrodinger equation using 
coupled-channel methods. We carried out direct least- 
squares fitting to the combined experimental data of this 
article and Refs. [14, 31, 36, 37, 82] to obtain a new 6- 
parameter model of the long-range interaction potential, 
which we designate M2012. 

The M2012 potential reproduces the experimental re- 
sults much better than the earlier M2004 potential [9], 
particularly at higher fields (above 250 G). It also pre- 
dicts j-wave and l-w&ve states, which were not included in 
the least-squares fits. The calculated positions of Z-wave 
bound states agree well with the experimental results re- 



ported in Ref. [15], for which the M2004 potential failed; 
this demonstrates the predictive power of the new model. 
The pole positions of the two broad s-wave Feshbach res- 
onances at high field are calculated to be 548.78 G and 
786.80 G for the -6(34)7s(6) and the -6(34)6s(6) states, 
respectively. 

The M2012 potential does have some remaining defi- 
ciencies. In particular, it is fitted only to results from 
ultracold collisions and the bound states that lies within 
10 MHz of the atomic threshold. It does not include re- 
sults from electronic spectroscopy on deeply bound levels 
of Cs2, or the near-dissociation levels observed in 2-color 
photoassociation spectroscopy [87], which are bound by 
5 to 80 GHz. It also does not satisfactorily reproduce 
the positions of i-wave states associated with excited hy- 
perfine thresholds. Resolving these remaining issues will 
require a simultaneous fit to all the experiments together, 
and is a topic for future work. 

Our new model allows us to make an accurate connec- 
tion between the experimentally controllable magnetic 
field strength B and the s-wave scattering length a over 
a wide range of fields. The scattering length is the essen- 
tial parameter in universal theories, and this connection 
was crucial in allowing us to interpret our measurements 
of three-body recombination in terms of universal Efimov 
physics in Ref. [18]. The present work is important not 
only for experiments on Cs, but also provides important 
information for ongoing and future experiments involving 
Cs mixtures, such as RbCs [88], LiCs [89, 90] and other 
interesting combinations. 
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